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1. Introduction 

The topological string theory has provided a lot of excitement in both physics and mathematics. 
It comes in two different versions, A- and B-models, based on the different twists of the worldsheet 
theory. Many techniques have been developed to compute their amplitudes. The techniques 
manifest different features of the topological string theory and are usually complementary to 
each other. For the present work, on the A-model side, the refined topological vertex [1] is the 
method of interest, and on the B-model side the remodeling approach [2] based on the topological 
recursion. We want to connect these two approaches to shed some light in the refinement of the 
mirror symmetry. 

The topological vertex formulation gives the all genus amplitude for the toric Calabi-Yau 
threefolds. This includes geometries which engineer the M = 2 SU{N) gauge theories in Ad 
in type IIA compactifications [3]. The genus zero amplitude of the topological string theory 
gives the prepotential of the gauge theory. Higher genus amplitudes compute the gravitational 
coupling to the gauge theory. The genus expansion of the topological string theory has been 
connected to the instanton calculus on the gauge theory side and allowed what is now called the 
refinement of the topological string theory. The refined topological vertex has been motivated 
to understand the A-model side of the extension. 

On the other hand, the connection between the topological string theory and matrix models 
has been under extensive exploration for a long time starting with [4, 5]. Both open and closed 
topological string amplitudes have been computed for various backgrounds using matrix models 
for both the topological A and B models. The integrability of topological strings on these 
backgrounds is one of the immediate consequences of this correspondence. More recently, the 
topological recursion is developed to solve the loop equations of a matrix model in a recursive 
manner [6]. The procedure starts with a spectral curve and basic ingredients associated with 
it. The spectral curve is obtained by the large N limit of the Schwinger-Dyson equations of the 
matrix model, then once the spectral curve is known, the topological recursion generates a set of 
meromorphic forms and complex numbers only based on the geometric properties of the spectral 
curve. Later the formalism is extended beyond the matrix model by choosing algebraic curves 
that are not coming from the large A^ limit [7]. Many remarkable properties are shown. 

The explicit relationship between the topological string theory on toric geometries and the 
topological recursion has appeared for the first time in [8]. It has been based on the observation 
that both the local B-model topological string theory amplitudes and the general solutions of 
matrix models can be related to chiral bosons living on a Riemann surface. The Riemann 
surface in question is the reduction of the Calabi-Yau threefold or the spectral curve, respectively. 
The mirror geometry of toric threefolds can be reduced to a Riemann surface as well. The 
meromorphic differentials are constructed on this mirror curve. This initial attempt is completed 
in [2]. This framework allows one to compute open and closed B-model amplitudes unambiguously 
at any point in the moduli space, but still perturbative in the string coupling. The integrals of 
the differentials correspond to the framed open amplitudes and the invariants are related to 
the closed amplitudes of the A-model, after taken the mirror map into account. However, this 
approach is restricted to the unrefined case. 

Our aim is to extend the remodeling to match the refined A-model computations. More 
precisely, we are after the refined spectral curve, i.e., the refined mirror curve, to perform the 
topological recursion to give the refined topological string amplitudes. We will first construct a 
matrix model for the refined topological vertex and the refined topological strip. Then we will 



find the spectral curves of these matrix models. The first step is straightforward considering that 
the refined topological vertex has been defined through plane partitions and a matrix model for 
plane partitions with arbitrary weights has been already constructed [9]. For the second step we 
will benefit from the results of [10, 11] where the spectral curve of the resulting chain matrix 
model has been found using saddle point approximation. We find explicit expression for the 
spectral curves as a function of /3 = — ei/e2. 

The closed refined strip computation has been connected to (chiral) Aat-i Toda theory (Tn- 
theory) three-point function with one of the three primary fields of a special type in [12] using 
the AGTW correspondence [13, 14]. Our construction for the refined strip implies a matrix 
model for the three-point function, hence, manifestly integrable structure for them. In [15] a 
analogous proposal for three-point functions of a 5d version Toda theory has been made based 
on the five-brane webs. The brane configurations are related to the toric spaces [16]. For T2, 
the corresponding geometry is the so-called closed vertex geometry and the chiral three-point 
function has been shown to agree with the 4:d limit of the refined topological vertex result. 
Higher rank theories involve more complicated geometries and it is not obvious how to take the 
4d limit. Our matrix model construction can be extended similar to [17] to capture all possible 
toric geometries. It will allow to connect matrix models to the proposal in [15]. 

The paper is organized as fallows: In section 2 we review the basic definition of the plane 
partitions and the description of plane partition in terms of non-intersecting paths. This descrip- 
tion is very crucial to write down the corresponding matrix model. In section 3 we overview the 
topological vertex formulation. In section 4 we remind the basic ingredients for the corresponding 
matrix model. We also make the connection to the refined topological vertex. In section 5 we 
extend our matrix model to cover the strip geometry in the refined case. In section 6 we focus on 
finding the spectral curve of the matrix models and the corresponding limiting shapes. In section 
7 we comment on the use of topological recursion. In section 8 we work out some examples and 
certain limit of them. In section 9 we finish with the conclusions. In the appendix we collect 
some useful formulae used in the derivations. 

2. Plane partitions 

In this section we want to review the basics about the 2d, also known as Young diagrams and 3d 
partitions, the so-called plane partitions. This section should serve more as fixing the notation 
rather than a detailed introduction. At the end of this section we want to remind how the plane 
partitions can be represented in terms of non-intersecting paths. In other words, they can be 
visualized as particles moving from one end to the other in discrete times and are only allowed 
to make jumps by ±1/2 to the next time slice such that no two particles can occupy the same 
position. 

2.1 Partitions 

Before defining the plane partition it is useful to review the 2d partition or the so-called Young 
diagrams. A sequence of positive integers of decreasing order jj, = {fii > //2 > Ms ^ • • •} defines a 
2d partition. By placing fii boxes on the i position the 2d partition can be pictorially visualized. 

/i = {6>3>2} — > I . 



It is also interesting to represent a 2d partition fi, rotated by 45 degrees, and define: 

hk = ^J^k- k + a, (2.1) 

where a is a fixed reference point, which we shall choose later on. hk is the horizontal coordinate 
of the end of row /i^ rotated by 45 degrees, see figure 1, i.e., the hk's correspond to the descending 
edges of the rotated partition. We have 



hi > h2 > hs > ■ ■ ■ > hjy > a - N. 
The profile of a partition fj, is the function r^{s) : M — t- M, defined by 



f2.2) 




h^h2 h^ 



Figure 1: The 2d partition rotated by 45 degrees, the horizontal coordinate of the descending edge is 
labelled by h^. The upper blue line is the profile of the partition. 

'^^J.is) = n Nl + T X] \s - {m - i + 1)\ + \s + i\ - \s - [fii - i)\ - \s + i - 1\, (2.3) 

i 

which is a piecewise affine function with slopes r' (s) = it 2, see figure 1. The slopes at integer 
s might be ill defined, but they are always well defined at half-integers. We define S-{^) to be 
the set of descending edges, i.e., the set of /i^'s: 



S-{n) = Is £Z\r'^{s- -) = --[ = |/ii,/i2,/i3,...| 



and we also define the set of ascending edges: 



5+(M) = |sGZ|r;,(s-^) = +^|=Z\5_ 



-(m)- 



(2.4) 



(2.5) 



2.2 Plane partitions 



A plane partition is the 3d generalization of Young diagrams. It is defined as an array of non- 
negative integers satisfying 

n,j > TTj+rj+s, r,s >0. (2.6) 



The plane partition can be pictorially represented by putting ttjj- boxes over each square 
labeled by {i,j). It can also be viewed as a crystal growing in the corner of a room: new boxes 
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Figure 2: An example for a plane partition, and its slices. 

can attach only in corners. A plane partition vr can also be represented by a set of 2d partitions 
vi^) — {^i(*) — ^2(5) ^ • • • }; labelled by s G Z. r]{s) is called the s**^ diagonal slice (see figure 
2), defined by 



■^ = {■■■, ^(-1), ??(0), r/(l), 7?(2), . . . }, where 7y(s) 



Vkis) = #{i I TTi,s+i >k} if s < 
Vk{s) = #{j I TTs+j,j >k} if s > ■ 

(2.7) 



2.3 Plane partitions and non intersecting paths 

Plane partitions are known to be expressed as a model of self-avoiding jumping particles. Consider 
a bounded plane partition, with three boundaries given by three partitions A, // and i^, located 



at respective distances Nx = N, N^ 



■Sinim -^^u 



from the origin, see figure 3. 



The plane partition can also be represented as the set of its diagonal slices rj{s), Smm < s < 
■Smaxi as defined in (2.7). Since vr has an height of A^, each slice r]{s) is a 2d partition with at 
most A^ non-empty rows. 

As we have seen above, a 2d partition 7]{s) = {?7i(s)}, with rows r/i(s) > r/2(s) > ••• > 



VNis) > 0, can be encoded by its 45 degrees rotated version, i.e., the set {/ij(s)}i=i,...,7v with 



hi{s) = r]i{s) -i + N + -\s\ 



(2i 



where we have chosen the origin (i.e. ?/Ar(s) = 0) to be on the Une h = 2\s\. hi{s) is the position 
of the i vertical edge of t/(s), along the s vertical line. 



The fact that vr has boundary A from above implies that 

hi{s) > h2{s) > ■■■> /ijv(s) > rx{s), 



(2.9) 



where rx{s) is the profile of the partition A bounding vr from the top, defined in (2.3). We obtain 
the 2d partitions by slicing a plane partition, this imposes the following condition among adjacent 
slices 



hi{s + l) = hi{s)± 



1 



(2.10) 



It is customary to represent the partition r]{s) by a set of N "particles" living on a Id 
line (coordinate /i), the position of the i^^ particle at time s is hi{s). Those particles follow 
a time evolution (the time is s, going from the left side to the right side of vr), such that 
at all times s the particles can never occupy the same position and can never cross (indeed 
hi{s) > /i2(s) > • • • > hN{s)), and move by jumps of zt^. See figure figure 4. 
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Figure 3: A plane partition with boundaries labeled by A, /i and i^. The red lines represent the condition 
(2.10). 

The trajectories of those A^ particles are N non-intersecting lines connecting the left bound- 
ary to the right boundary. The intersection of the lines with the (vertical) time lines are the 
location hi(s) of the particles at each integer time s = Sminj • • • ; Smax- Each of the A^ lines will 
represent the world-line of a particle discretely moving in one dimension, see figure 4. 
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Figure 4: For each plane partition there is a set of trajectories of particles moving in Id from the left 
boundary to the right one, subject to the boundary conditions. 



A plane partition configuration, is thus encoded by the data of the positions hi{s) of A^ 
particles (i = 1, . . . , A^) at times s with Smm < s < Smax, whose initial positions /li(smm) and 
final positions /li(smax) are fixed: 



vr = {/li(s)}j=i,...,Ar 



,<S<£ 



such that 



hi{s) > /i2(s) > 
hi{s + 1) = hi{s) 

'ijl.'Smaxj — l^i 



> hN{s) > rx{s), 



2' 



i + N + 
■i + N + 



2 ■ 



(2.11) 



3. The topological vertex formalism 

The topological string theory offers a simplified framework to understand the physical string 
theory. It shares many important phenomena with the string theory but in a computational 
more tractable manner. The topological string amplitudes appear in the Ad effective action in 
string theory compactifications. They are also generating functions for BPS states in various 
supersymmetric theories. 

Supersymmetric gauge theories can be constructed within the framework of string theory 
in different ways. One of the most powerful constructions is the so-called geometric engineering 
[3]. The gauge theory data such as the gauge group, the matter content and the amount of 
supersymmetry is encoded in the geometry of the space X that is used to compactify type IIA 
string theory from Wd down to Ad or M-theory from lid to 5d. The geometric engineering gives 
a recipe to determine this space X based on the gauge theory data. One of the most important 
advances in string theory is the realization that the genus zero amplitude of the topological string 
theory on X computes the prepotential of the corresponding gauge theory. Higher amplitudes 
compute the gravitational corrections. More precisely, in type IIA compactification the effective 



action in the 4d transverse space includes F-terms describing the couphng of the self and anti-self 
dual graviphoton fields to self and anti-self dual contractions of the Riemann tensor 

f d^x Fg{co){RlFl'-^ + RlF^'-^), (3.1) 

where Fg{oj) is the genus g topological string amplitude. It depends on the Kahler moduli 
Tp = jo 00 of the internal space. They are generating functions for genus g Gromov-Witten 
invariants and are computed by performing integral over the moduli space of Riemann surfaces. 
Turning on only the self-dual graviphoton field we get corrections to the i?^ term 

oo \ 

yFl^-^F,{u:)\R\, (3.2) 






where F+ plays the role of the topological string coupling constant. The partition function is 
defined 



oo 



Z(^,ff,) = e^(-'^^) = exp Y^al'-^^si^) ■ (3-3) 

The target space point of view offers a more physical interpretation of the topological string 
partition function [18, 19]. If we compactify type IIA on X, in the remaining 4d space we 
will have BPS particles which come from D2 branes wrapping two cycles /3 G H2{X,Z) in X. 
The mass of each particle is given by the area of the cycle T^ the brane is wrapping. There 
is a caveat due to the bound DO branes to D2 branes. Let us lift to the M theory, then a D2 
brane becomes an M2 branes wrapping the corresponding two cycle. An M2 brane comes with 
momenta along the circle direction. We can interpret these Kaluza-Klein modes in the M-theory 
with the DO branes. The mass of the M2 brane with the momentum n is shifted from T^ to 
Tp + 2T:in/ gs- In addition, it carries quantum numbers under the little group of massive particles 
in 5(i, {jl,Jr) G SU{2)l x SU{2)fi = 5*0(4). Let us denote the degeneracy of the particles of 
given quantum numbers by N^'"''"^ . The free energy can be written as contributions coming 
from all such particles: 

in terms of N o^ = ^ • (— l)^-'«(2j/j -|- 1)N£^'''"\ not in terms of the individual degeneracies. 
This is because the individual degeneracies may change as we vary the complex structure due to 
pairing of different Jr into multiplets. There exists an extension of the topological string theory 
due to the advances in the microscopic derivation of the prepotential of the supersymmetric 
gauge theories [20, 21, 22]. This extension is equivalent turning on both the self and anti-self 
dual graviphoton fields [23]. The free energy turns out to be written in terms of the individual 
degeneracies 



F(uj:t, 



Z^ 2-^2-^ „(fn/2 _ f-n/2\(^,n/2 _ „-n/2\ 



riTf, 



(_l)2iL+2jfi^&J«) ('(tg)-«iL + . . . + (tg)"J'i) ((p""^« + . . . + CqPA 

(3.5) 



where q = e + and t = e . 

Our main interest is AA = 2 supersymmetric SU{N) gauge theories with or without matter. 
This type of theories are engineered using the toric geometries. The toric geometries are encoded 
in Id tri-valent graphs which show the degeneration loci of the toric action on them in terms of 
the slopes of the lines in the toric diagram. The topological vertex formalism offers an all genus 
computation of the topological string partition function for the toric Calabi-Yau threefolds. 





Figure 5: The C^ patches. The patches can be glued to give rise to general three-fold, e.g., 0{—l) -^ P^. 

The idea behind the topological vertex [24, 25] is to divide the geometry in C^ patches, 
as shown in figure 5, and compute the open amplitudes for the individual patches. Different 
type patches arise in cutting the geometry but they all are related to each other by symmetry 
transformations of the theory. This allows to choose one generic patch and define the topological 
vertex for it. From the point of view of string theory this cutting is realized by placing Lagrangian 
branes/anti-branes. The topological vertex C\^y{q) is defined through the open amplitude of C^ 



^(^1,1^2,^3) = Y. Cxf,u{q)tTxVliT^V2iTuV^, 



(3.6) 



X,fj,,iy 



as a function of the parameter q = e~^' with gs = (-F+) being the topological string coupling 
constant. T^'s are the holonomies supported on the branes, and the sums are over irreducible 
representation of U{oo). A suitable gluing algorithm joins the amplitudes of the individual 
patches into the amplitude of the whole geometry. 

A combinatorial interpretation the topological vertex is given in terms of plane partitions, it 
is often called the "crystal model" [26, 27]. It is identical, up to the MacMahon function, to the 
generating function of plane partitions subject to boundary conditions dependent on the repre- 
sentations labeling its legs, see figure 6. The (formal) counting parameter on the combinatorics 
side is identified with the exponential of the topological string coupling constant. The generating 
function is computed by dividing the plane partition into slices, see (2.7). 

The refined topological vertex is originally constructed modifying this combinatorial interpre- 
tation of the usual topological vertex to reproduce the equivariant instanton part of the partition 
function. For the refined topological vertex, the direction along the slices of the plane partition 
plays a significant role unlike the usual one. The representation along this so-called preferred 





(a) 



(b) 



Figure 6: (a) 7r,(A,^,i/) for A = (6,4,3,1,1),// = (5,4, 3, 2, 2), z^ = (4,3,2,1). (b) An example of 
7r(A,/i, i^). The topological vertex is a generating function counting configurations 7r(A,/i, i^) for given 
A, /i, v^ with a weight depending on slices of the plane partition. 




■K — ^ 



n^ gi"^"^l = g^-1 l-(-*-^)l i^T=i I-' 



(-'^j+i-i) 



g = blue (solid line), t = red (dashed line) 
z/=(4,3,l) 



Figure 7: The plane partition seen from above. The grey zone is the partition A. Slices at times s E S- (A) 
are weighted by il''(*)l, and slices at times s £ S+{X) are weighted by gl''('*)l. 

direction determines the weights of the slices, i.e., which slices are counted with q or t. This is 
depicted in figure 7. 

Consider a plane partition vr = {^(s)}, and its slices r](s). The refined topological vertex is: 



vr, 57r=(A,/x,i/) se5-(A) seS+{X) 



(3.7) 



It is well known that if we have no boundary X = fi = v = 9, and t = q, this is the 
Mac-Mahon function 



^0,0,0(9,9) = n^^"^ 



k\-k 



(3.8) 



k=l 
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3.1 Weight of a plane partition 

The bounded refined topological vertex, is the generating function counting plane partitions 

^ = E ^(^) (3.9) 

with the weight (3.7): 

W{tt)= Yl t^.-i^'W Yl g^^i^^^"^. (3.10) 

seS-{A) ses+(A) 

Notice that the restriction /lAr(s) > rx{s) can also be encoded by assigning a vanishing weight to 
plane partitions which do not satisfy that condition, and thus we write 

mvr)= n n^(^^(^)'^)' (3.11) 

where 

r t'', if s £ S-{X) and h - rxi-s) G N 

w{h, s) = < g'*, if s G S+{X) and h - rx{s) G N (3.12) 

(o, if h-rx{s)(^n. 

More generally, one could generalize this problem, and require that the possible paths are 
restricted to a certain domain T> in the (h, s) plane, of an arbitrary shape (not necessarily given by 
a partition A), or even more generally, the paths could be weighted by an arbitrary "landscape" 
potential in the (/i, s) plane, i.e., one would like to give an arbitrary local weight w{h, s), so that 
the weight of the path is the product (3.11) and where the weight function w{h, s) can be chosen 
at our convenience, for example a forbidden position is such that w{h, s) = 0. 

In our case, for the refined topological vertex, non-forbidden positions will be chosen with 
weight w{h,s) = t^ or q^ depending on wether s is in the ascending or descending side of the 
upper partition A, i.e., if 'w{h, s) / 

, N fi^, when sgS'_(A) . . 

^(^'^) = l,^ when sGS^iX)- ^^'^^^ 

Notice that there is no need to assign values to ■w{h,s) at positions {h,s) which can never 
be visited by the paths (notion of shadow of a domain in [9] ) , and therefore we may also choose 
w{h, s) = t^ or q^ in a forbidden region which will never be visited. 

This implies that we need to choose a non-trivial weight, only for s being an outer corner of 
A. We write: 

5_(A)= U]ui^i,Vi] , S+iX)= U]vi,Ui], (3.14) 

1=1 1=1 

i.e., Vi are the inner corners of A and Ui are the outer corners of A, see figure 8. Eventually, for 

the refined vertex we shall use: 

' t'^, when Ui-i < s < Vi 

^<l' (l-Elt!r''^'"' -a::!;S5)" ) /(^-^a(^-))> when . = .. 
where f{h) can be any function of h whose value is 1 when /i G Z. We shall write 

w{h,s)=e-^'^^'\ (3.16) 

and we call Vs{h) the "potential" at time s. 



11 



4. The matrix model for the refined topological vertex 

In this section, first, we want to review the general construction of the matrix model for plane 
partitions with different weights along each 2d slice when we slice it. Later, the resulting matrix 
model is restricted based on the choices made for the refined topological vertex. It turns out 
that the potentials along most of the slices, more precisely, every slice but the first, last and 
along the outer corners, are linear. The simplicity of the potentials allows us to group slices into 
building blocks between the slices along the outer corners which we can glue afterwards. Within 
the individual blocks all matrix integrals but one can be performed leading to a significant 
simplification of the model. 

4.1 The matrix model for the plane partitions 

In [9], it has been proved that the sum over plane partitions with arbitrary weights w{h,s) can 
be rewritten as a matrix integral. Namely, theorem 3.1 in [9], says that the generating function 
of plane partitions is proportional to the following multi-matrix integral 

2 ^max 

^^ g-Try.(M.) (4_-L^ 

s' = s ■ -I-- 



/» -^max J- /• ^^^^.- ^ ''max 

/ 11 dMs jT dRs' n 



2 *m!ix 2 



-TrC/,,(fi.') TT ,TrR,,{M^,^i-M^,_i) 

a — -^min i^ 9 



X n e-^^^^'(^^'^ n e^''"='^"^»'+r^''^='-i\ (4.2) 



2 



where the potentials Vs{h) encode the weights w{h, s) = e ^''^^' and the potential Us'{r) encodes 
the jumps by ±1/2, i.e., e~^''^'^' is the Fourrier-transform of 5{h — \) + S{h + |), namely it is 

g-f/,/ (r) ^ g-mr ^ ^inr ^ ^ COS TTr . (4.3) 

The boundary conditions at the slices — Smin and Smax are encoded in the following fixed diagonal 
matrixes which are not integrated over 

^s^in = diag(/li(Smin), /l2(Smin), • • • , /iAf(Smin)), (4.4) 

^Smax = diag(/li(Smax), /i2(.Smax), • • • , /iAr(Smax))- (4.5) 

The proof of that formula is given in [9], but let us just sketch it here. It uses the fact that the 
generating function of non-intersecting paths is a determinant (the Lindstroem-Gessel-Viennot 
determinant LGV [28]) 

Z= det Wii (4.6) 

where Wij is the weight of an individual path going from /li(sinin) to /ij(sniax)) ^-e.. 



"Jj 



/OO ^max 
JJ w{Xs, s) dXs 5{Xs^i^ - hi{s-^i^)) 6{Xs^^^ - /ij(Smax)) 
-OO „_„ . 

^ -' ' 1. ., 1 



(4.7) 



X 

S — Sniin 



JJ ( 6{Xs+l - Xs - -) + 5{Xs+l -Xs + 
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The (^-functions appearing in the weights can be Fourier transformed 

S (xs+i -Xs-l) + 6 (xs+i -Xs+^^= j" dr e2-K-.+i-.) [g— + g-^] _ (^^^ 

Inserting the Fourier transform into the LGV determinant yields (after symmetrization over the 
integration variables) the determinants det e^*'^''*''^, which we recognize as the Harish-Chandra 
Itzykson-Zuber integral [29, 30]: 

det (e^'^ ''^''A = A{h) A{r) [ dt/ e^" ^ '^^^^^ (4.9) 



where h = diag(/ii, . . . ,/iAf), r = diag(ri, ...,rN) and A(/i) = l\iyjihi-hj), A(r) = nj>7( 



n-ri 



are the Vandermonde determinants. This allows to rewrite the sum over self-avoiding paths as 
a matrix integral. The eigenvalues of the matrices are the /ij's and the r^'s for each time s. 
More precisely, Mg is a matrix whose eigenvalues are the hi{s), and Rg,i is a matrix whose 
eigenvalues are the Fourier variables corresponding to the delta functions encoding the jumps 
from s —7- s + 1. The angular parts of those matrices are provided by the U G U{N) in the 
Harish-Chandra Itzykson-Zuber integrals (4.9). 

We shall not enter further details of the derivation of formula (4.1) here, since we are only 
interested in an application of this general formula, and thus we refer the interested reader to [9] 
for the full derivation. Another presentation of the derivation is also explained in [31]. 

4.2 Integrating out the intermediate matrices 

The matrix model (4.1) is a multi-matrix integral, known as the "chain of matrices" (see [32]). 
This is an integrable, exactly solvable matrix model (see [32, 33]). However, in our case, it can 
be further simplified. Indeed, for most values of s {i.e., s ^ Ui), the potentials Vs{Ms) are of the 
form 

Vs{Ms) = -Mslnq or ^^(M^) = -M^ Int, (4.10) 

i.e., the potentials are linear in Mg. Most of the integrals over Mg can be performed explicitly, 
and our multi-matrix model reduces to a matrix model with fewer matrices. The only Mg which 
can not be integrated easily, are those corresponding to values of s where the potential Vg{Mg) 
is not linear, that is when s = Uj is at an outer corner of the upper partition A. 

Therefore, for any given partition A we will group slices from one outer corner to the next 
one into a basic building block. Let us call the block from the (z — 1)* outer slice to the i as 

„ Ui — l „ Mj-l 

^.-l,(M„^_,,i?„^+i,M„J^ / , n ^^W r n n '^^^' 

s'^Vi+l/2 
Tri?,/(M^/+i/2-M,,_i/2)^ (4_11) 

The appropriate weight to compute the refined partition function is encoded in the potentials 

VgiMg): 

IntMg, for Ui-i < s < Vi 



Ur-l 


Mi-1/2 


Ui-1/2 


^ ^-TrVs{Ms) 


g-Trt/,,(i?,,) 


e 


S=Mi_l + l s' 


=-Uj_i + l/2 s' 


=«,_l+l/2 



^ ^ InqMg, for Vi < s < Uj. 
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V2 U2 V3 




Figure 8: The choice of weights. Red slices (i.e. ascending edges of A) have weight t, and blue slices have 
weight q. 



The integrals along the slices between the outer corners can be performed giving rise to matrix 
5-functions 



„ Ui — l Ui — ] 

/ n ^^^« n 



Ui-l/2 
-TrVsiMs) TT gTrR^,(Af^,_^i/2-Af^,_i/2) 

Mi_l + 1 S=Ui_l + l s'=Mi_i + l/2 

Vi-l/2 Ui-3/2 

n KRs'-Rs'+i+^^t) JJ 5(iiy-i?,.+i+lnQ)e'^'^^^"^^">-i/^"^"-i^"-i+^/2^ 

s'=«,_i+l/2 s'=Vi+l/2 

(4.12) 
After performing each Rgi integration but -R„.+i/2 in the interval ]ui-i,Ui[ one obtains 

Z^ (M R M ) = gTr-R(Mui-M„-_-^) (t),-«,_i)lntTrM„._j+(wi-i;i-l)lngTrAf„; 



Vi—Ui-i 

n 



-TrC/(/?-fclnt) 



fc=l fc=l 



n 



-Tr;7(iJ+(A:-l)lnQ) 



(4.13) 



The last two products can alternatively be written as 



Vi-Ui-1 



\[ e-^rU{R-kXnt) ^ ^{v. 



k=l 

Ui—Vi 



n 

fc=i 



-TrU(R+(k-l)lnq) ^ -( 



oc e 



-ni_i)Tr R/2 dct ( 1 6 \t) ^ 

det(-t^'-«i-i+ie-^;t)^ 

(«i-i),)Tr fl/2 det (-e ]q)^ 

det (-g("'-^')e^;g)^ 
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up to factors independent on R. (a; q)oo is the Pochammer symbol and is defined as 

oo 

ia;qU^ll{l-aq''). (4.14) 

fc=0 

The individual blocks thus take the form 

Tri?„^+i/2(Mu^-M„-_-^) (i;i-Ui_i)lntTrM„;_j+{«i-i)i-l)lnqTrM„ (2Di-«i-u,_i)TriJ„^+i/2/2 



--i— 1,; 



(4.15) 

det ( -t e~'^"«+i/2 ; t) det f-e^^.+i/^ ; q) 

(4.16) 



det ( -r«-««-i+i e"^"«+i/2;t ) det (-^K-^O e^"«+i/2; 



The whole integral is then given in terms of these blocks, hence, we end up with the following 
chain matrix model for the refined topological vertex 



Z= / lldMu, / di?,^+i/2e-^^^-'("'"-)^^-i,(M„,_„i?,^+i,M,J. (4.17) 



Let us remind that the representations of the refined topological vertex along the unpreferred 
directions are imposed by the fixed matrices at the first and last slices. The representation can 
be alternatively given using the inner and outer coordinates which appear in the above matrix 
integrals. 

5. The matrix model for the refined strip geometry 

In this section we want to extend the matrix model for the refined vertex to more general cases and 
propose a crystal model for the refined strip geometry [34, 35]. Two pieces of these geometries 
can be glued along the external legs to engineer interesting SU{N) geometries with extended 
supersymmetry, including 4d superconformal theories. 

Our model will resemble the bubbling picture of [36]. We will excise a region such as the one 
in figure 9 whose sizes are related to the Kahler parameters in the strip geometry. We grow the 
crystal in the remaining unbounded region. This geometry requires to redefine the weights on 
the slices along the outer corners to get a combinatorial model, each such slice is counted with t 
instead of q like in the case of the refined topological vertex. We want to work out the example 
shown in figure 9 in detail. From our approach it will be clear that there is no obstruction in 
generalizing this example to a longer strip. 

We will compute the generating function for the crystal using the transfer matrix approach. 
The following identities are crucial for our computation. The vertex operators r-i-(l) and the 
Hamiltonian Lq satisfy 

T.{IM = Y,W), T+{l)\fi) = Y,H, and g^1/x) = gl^'lA^). (5.1) 

We will compute the generating function in the following way: we will assign an arbitrary 
Young diagram to the slices along the inner and outer edges and find the contributions with these 
'boundary' conditions. In other words we are inserting identity operators at certain places in the 
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Figure 9: (a) the definition of slices, (b) tfie region where the crystal grows, (c) the refined strip geometry 

generating function to ease the computation. Later we sum over all possible Young diagrams 
along these slices and match these sums to the summations over Young diagram we obtain from 
the refined topological vertex gluing. It will be obvious that the continuation of this procedure 
iteratively and straightforward. The generating function for the example given in the figure 9 is 
given by 



oo ni — 1 mi—m2 — l 

z= E (0in'?^"r+(i)|A)(A|t^°r_(i)n*^"r_(i)|^)(<T|t^or+(i) J] 9^°r+(i)|^) 

X,a-,u,iJ,,r] i=l i=l i=l 



"2— '11 — 1 



(Ht^^r^i) n t''"r_(i)|^)(M|t^''r+(i)|r?)(7?|t^onr-(i)*^°|0) 



(5.2) 



j=i 



?=i 



Note that the skew Schur functions can be written as 



Sx/f,{xi,X2,...,Xn,0,0,...) 



E 



JAh|A(l)| |A{1)|-|A(2)| |A(n-l)|-|;.| 



(5.3) 



A!-A(l))-..>A(n-l)^^ 



Using the above identity and others collected in the appendix after some algebraic manipu- 
lations the generating functions takes the following form 

Z = M(t q)^ y^ (_Qm2y«i| /•_^n2-niya2| /■_„mi-m2\|a3l /_^»ii-)|a4| ^|ai|+|a2|+|a3l+|a4| 
cr,iJ,,T,{ai} 

X g-l/^l-l'^lt-l^l Sa, (1, t, . . OSa'^ll' 9' • • ^^^^/.(l) t, . . .)Sa^/^il, q, . . ■)Sa,.,/r{'i^, t, . . .) 

X SayM,q, . . .)s^tj^{l,t, . . .)saA'^,q, . . .). (5.4) 

This expression looks similar to the sums we obtain when we glue refined vertices with 
the choice of the preferred direction along the external legs. Let us briefly review the refined 
topological vertex computation. The partition function is given as 
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Z= Y. (-Qi)'^'(-Q2)l"l(-Q3)l^l(-Q4)l"IC^00(i,^)C^*.0(g,i)CA.*0(i,^)CAt,0(g,t)Q.*0(t,g) 
X s^(l, t, . . .) V/m (1' 9' • • ■)^y'/-ni (1' *' • • ■)^y/m (1' 9, • • O'SA/r^a (!'*>•• O^AV^ (^i 9> • • •) (5-5) 

X Srt/r^,^{l,t,...)Sr{l,q,...) 

The refined vertex result matches the crystal model under the identification of the Kahler 
parameters 



Qlyf = g'"^ Q2^j = e'-''\ Q3^^ = g'"^-'"^ ^4^^ = ^^ (5.6) 

up to MacMachon fucntions 

^vertex ^ J^yt^Q) ^crystal- V^-') 

Having found a crystal model for the refined strip geometry it is straight forward to derive 
the corresponding matrix model. As mentioned the counting scheme differs from the one of 
the refined topological vertex at the outer corner. Luckily this difference can be very easily 
accommodated in our previous computation. We integrated out the intermediate slices between 
two outer corners and defined the individual blocks Zi-i^i. Since there are no modifications for 
these intermediate slices the blocks remain unchanged. However, when we are computing the 
partition function we need to integrate with respect to these modified slices. Therefore we only 
need to change the potential along these slices 

!■ M 

Z= WdMuJ di?,^+i/2e-T'-^"^(^"«)^,_M(M„,_,,i?„^+i,M„J, (5.8) 

with K, (M„J = - In tM„,. 

6. Spectral curve and limitting shape 

Our next aim is to find the spectral curve of our matrix model. It was derived in [10, 11] that 
the spectral curve of a chain of matrices is obtained as the "saddle point" of the corresponding 
matrix model. Before we move to the details, we want to point out a caveat of the saddle point 
approach. Unless the potential is polynomial it is possible that there are infinitely many saddle 
points of the potential. The spectral curve is the one which gives the global minimum, hence, 
for such potentials it is almost impossible to check whether this condition is satisfied. The lack 
of this explicit verification of the global minimum is not a big hindrance to obtain the spectral 
curve. However, when (3 is rational, all those conditions can be realized with algebraic functions, 
and the natural guess, is to look for the smallest possible degree algebraic curve which satisfies 
all the conditions. 

Let us first review the basic strategy to find the spectral curve using saddle point approach. 
Later we motivate and write down our algebraic ansatz. Finally, we comment on limiting shapes 
of the crystals. 
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6.1 Rules for finding tiie spectral curve 

According to [10, 11], the spectral curve is described as follows: there must exist a Riemann 
surface C (to be determined later from consistency conditions), and some analytical functions 
C —7- CP^ (one per matrix in the chain), which we note: 

Z>-^Xsiz), Smin < S < Smax (6.1) 

Z^rs>{z), Smm + - < s' < Smax- - (6.2) 

such that the saddle point conditions read as 

rs~i/2{z) - rs+i/2{z) = V^{xs{z)) (6.3) 

Xs'+i/2iz)-Xs,^i/2iz) = U!.,{rs'iz)) = -- ^^^^(^) (6.4) 

While deriving our matrix model we divided the times into blocks based on the representation 
along the preferred direction. Similarly, we will treat these regions separately: 

• For slices in the region Uj_i < s < Vi the potential is 

Vs{x) = -x\nt. (6.5) 

This allows us to define a function Ri{z) = r^.i{z) such that 

Vs, Ui-i<s<Vi, rsj^ii2{z) = Ri{z) + {s -Vi)\nt, (6.6) 

and thus 

xs+i{z) - xs{z) = -- ^,,,^^^^(,)^^ - (6.7) 

Let us define: 

Yi{z) = e^'^^\ and X,(z) = t"^-(^). (6.8) 

Summing over some intermediate slices 



Vj-S 



1 Z^YAz)-t^ s-v, ^ t 



ffc 



X.M) - xsiz) = -- g ^^^y^, =^ + J:^ WfT¥- ^'-'^ 



In the limit In t — )■ this gives 

Xsiz) t-^+t-^Yi{z) 



VA—S S — Va 



(l + 0(lnt)) (6.10) 



X,^{z) l + Y,{z) 

Similarly, for the region Vi < s < Ui we have 

Vs{x) = —X Inq. (6.11) 



this implies that 
and thus 



rs+i/2iz) = Riiz) + {s - Vi)lnq (6.12) 

I qS-Vi ^Riiz) _ 2 

Xs+l(z) — Xs(z) = f7^^ . (6.13) 

+^ ^ ^ ' 2 qs-'v^QR^{z) _^i ^ '' 
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In the same way we have 

Let us define again 

X,(z) = t-^^("). (6.15) 

In the limit Int — t- this gives 



where we have introduced 



(l + 0(lnt)) (6.16) 



f3 = lnt/lnq (6.17) 

that is assumed to be 0(1) when Int — )■ 0. 
• At times s = Ui, we have 

r„^_i(^)-r„^+i(^) = <(x„,(z)). (6.18) 

Remember that the potential Vu^ix) is in this case not well defined. The only requirement 
is that it is worth Vui{x) = —x Inq if a; is an allowed position, and Vu^ix) has a logarithmic 
singularity at all non-allowed positions. This means that within the allowed domain, we can 
replace Vui{x) by —x Inq, and at the boundary of the allowed domain, it has a singularity of the 
form ln(X„^(z) - T^K)) + O(lnt). This implies that 

yi{z) 

Yi{z) = 0, ifX„^(z) = f^(«») (6.20) 

Y,iz) = oo, if Xu,_, (z) = r^("-i) (6.21) 

6.2 Ansatz: an algebraic spectral curve 

We see that to leading order in Int, we need to find functions Yi{z) which are analytical on 
a certain Riemann surface C, which have a zero at Xui{z) = f^^"') and a pole at Ar„._-^(2;) = 
f\(ui^i) ^ Moreover, we need to find functions X^.{z) analytical on C which vanish at Yi{z) = —1 
with exponent 1 and /3, and which diverge at 1^ = —fi-'^i-'^ with exponent 1 and diverge at 
Yi = —q^i~'^i with exponent /3, and such that 



V^—S S — Vj 



Xs{z) = X^^{Z) *^]^*y^^'^''-' , Ui-l<S<Vi 

Xs{z) = X,Xz) (l^^^^^I^\ , v,<s<ui (6.22) 



^i+l(^) — 4-Vi + l-Ui „Ui-V, 



Yiiz) 



^v,+i-u, qUi-v, ^ S = ui. 



All those requirements are of algebraic nature, and it is natural to look for an algebraic curve 
C and algebraic functions of the minimal possible degrees satisfying all the requirements. 
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It is rather easy to find this minimal degree algebraic curve, as one can see in the examples 
of section 8 below. As we mentioned, there is no proof that this minimal degree algebraic curve, 
is indeed the spectral curve of the matrix model, i.e., that it is indeed the leading saddle point 
corresponding to our matrix integral, and that there is no other saddle point with lower energy. 
However, in many examples, this choice looks very reasonable. For the unrefined case t = q, 
this minimal degree algebraic curve coincides with the one found by Kenyon-Okounkov-Sheffield 
[40]. See an explicit example in section 8 below. 

From the point of view of string theory, that minimal degree algebraic curve, in the unrefined 
case, is the mirror curve of the geometry under consideration. We thus propose that, in the refined 
case, our minimal degree spectral curve is the mirror of our refined geometry. 

6.3 Limiting shape 

In [10], the functions Xs{z) and rs{z) have been introduced as the limit resolvents of matrices 
Ms, namely, if one eliminates z between Xs{z) and rs(z), we have almost by definition 

Tr . . ,^ ) ~^+i(^), (6.23) 

Xs{z)-Ms/ ^2 

where the upperscript {.)^^' means the planar part of the expectation value, i.e., the leading Int 
term. On the other hand, the resolvent is the Stieljes transform of the density ps of eigenvalues 
of Ms which means that the density of eigenvalues of Ms is the discontinuity of the resolvent 
across its cut 

Psixsiz)) = ^ (r,+ i (z_) - r,+ i (z+)) (6.24) 

where z_ and z^ mean the value of z on the right and on the left of the cut. 

By construction, the eigenvalues of Ms are the hi{s) defined in (2.8), and they represent 
the heights of the slice s of the plane partition, or in other words, hi{s) is the position of the 
i particle at time s. Equivalently, it is the position of the i non-intersecting line at time 
s, hence, the position of the i^^ vertical edge of the plane partition vr at time s. Ps{h) is the 
density of particles at time s, i.e., the density of non-intersecting lines at time s, corresponding 
to the density of vertical edges at time s. See figure 10. Therefore, from the knowledge of the 
functions Xs, rs, we can recover the limiting shape of our plane partition. This has been done 
in a general setting in [9], and corresponds (in the unrefined case) to the limit shape found by 
[37, 38, 39, 40, 41]. 

6.3.1 Reminder unrefined case 

In [37, 38, 39, 40, 41], it has been found that the limit shape of the plane partition is such 
that there is a "frozen" region (where the density of particles = vertical edges is constant) and 
a "liquid" region. The curve separating the frozen region from the liquid region is called the 
"arctic circle", see figure 11 and figure 12. 

Furthermore, in [40], it has been shown that the arctic circle is a real algebraic curve, 
inscribed in the domain of the plane partition (the polygon bounded by the profile rx{s) on the 
top, and by the lines s = Smm and s = Smax on the sides, and by the lines h = N — Smin + s/2 and 
h = N + Smax — s/2 at the bottom), tangent to all sides, and with the minimal possible degree 
given by 

degree = ^ (5 + 2#{n,} + #{v^}) (6.25) 
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density 



\ih=xJz) 



Figure 10: From the knowledge of the functions Xg, rg (i.e. the spectral curve) one can recover the 
density Ps{h) of eigenvalues of the matrix Mg, and from there the limiting shape of the plane partition 
can be obtained. The arctic circle is obtained by Legendre transform of the limiting shape. 



min 




max 



N 



Figure 11: The arctic circle is tangent to all sides of the domain, and has cusps at outer corners of A. 

For example, if A = is the empty partition, the profile r\{s) = \s\/2 has only one inner 
corner vi = 0, and no outer corners, i.e., ij^{ui} = and #{^4} = 1. In this case, the arctic circle 
is an algebraic curve of degree 2, that is the unique ellipse inscribed in the hexagon. See figure 
11 for the typical shape of the arctic circle. 
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Figure 12: The statistical behavior of the surface in the liquid region is governed by the Burger's equation. 



In general, the arctic circle found in [40] is called a Harnack curve and has many nice 
properties, see [40] for details. In addition, it has been shown that, within the liquid region, the 
small Inf statistical behavior of the surface of the plane partition is governed by the Burger's 
equation. See figure 12. We believe that the loop equations of the corresponding matrix model, 
are equivalent to the Burger's equation, although this has never been carefully proven. 

6.3.2 Some general remarks about the limit shape 

The general properties of the arctic circle curve have been studied for any weight w{h, s) in the 
domain in [9]. In particular, it has been shown that the arctic circle is always tangent to the 
boundaries of the domain, and has vertical tangents at times s = Ui, i.e., at outer corners of 
the boundary. For any weight w{h, s) the limit shape is always continuous and has a continuous 
tangent at everypoint except at cusps. In particular for our case, this means that the arctic circle 
is continuous and has a continuous tangent at inner corners s = Vi. 

It has been also shown in [9] that the arctic circle curve and the spectral curve {xs{z),rs{z)) 
are related by Legendre transform, namely, let h = x{s) be the equation of the arctic circle in 
the (/i, s) plane. Remember that the spectral curve is 



Xs{z)=t-^^^'^=X,^{z) 



+ t 



2 



n^) 



l + Y,{z) 

of the general form 

Xs{z)=t^ {A{z) + t'B{z)), 

i.e. t^''^Xs{z) is linear in t^. The value at the extremum given by dXs{z)/dz = is thus 



t 



-x(z) ^ ^^ (extremum {A{z) + t" B{z)) 
i.e., it is a Legendre transform in the variable i*. 



t 2 Legendre(A o i? ^)(t 



(6.26) 



(6.27) 



(6.28) 
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The reverse Laplace transform can also be explicitly performed as well, let t^ = t^^''^ Pit^) 
be the equation of the arctic circle in the (/i, s) plane. We can recover the spectral curve 
{xs{z),rs{z)) by choosing the spectral variable z to be (locally) the argument of F, namely 

= t-^/2 (F(z)-t'F'(z)) 

^ ^ ' ^ " (6.29) 

Note that equations (6.20) and (6.21), are precisely the conditions that the limit shape is 
tangent to the boundary of the domain. 

7. Topological recursion and mirror symmetry, "remodeling the refined B- 
model" 

The chain of matrices has been completely solved in the following sense: 

• We can write loop equations, i.e., identities among correlation functions are obtained by 
integration by parts. Loop equations for general chain of matrices have been written in [10]. 

• If we assume that all correlation functions have an expansion in some small parameter 
Qs, such that the connected part of the n— point function behaves as gs^~ ", then the general 
solution of loop equations for the chain of matrices has been derived in [11]. It has been found 
in [11] that the solution of loop equations is given by the topological recursion defined in [7], 
applied to the spectral curve we have just described in the previous section. 

The axiomatic definition of the topological recursion has been introduced in [7]. It is a 
recursion starting from the data of a Riemann surface C and two analytical functions x : C — )■ 
CP^, y : C ^ CP^. It associates some sequences of symmetric meromorphic differential forms 
WJf^ on ®"r*(C) =T*{C) (g)---(g) T*{C), and some complex numbers Fg{S) = WJ^^^ G C to a 
spectral curve S = {C,x{z),y{z)}. The recursion goes on — x = 2g + n — 2, i.e. Wn is computed 
in terms of Wl^, with 2g' + n' — 2<2g + n — 2 i.e. — x' < —X- We refer to [7] for details. 

The numbers Fg{S) are sometimes called the symplectic invariants of S. The main result of 
[11] is that the formal small gs = Int expansion of the logarithm of the partition function of a 
chain of matrices, is given by the symplectic invariants of its spectral curve 

oo 

InZ = Y,{^ntf'~^ Fg{{lnXs,lnY}), (7.1) 

9=0 

where Fg is the g^^ symplectic invariant of the spectral curve {C, In Xs{z), In Y{z)} (and it is 
noticed in [11] that symplectic invariance implies that Fg{{C , In Xs,\nY}) is independent of s, as 
expected). In other words, knowing the spectral curve, one can recover the full Int power series 
expansion of the partition function, by applying the topological recursion to the spectral curve. 

We shall not enter detailed computation of symplectic invariants of the spectral curve we 
have just found, but we just mention that this statement is the refined version of the "re- 
modeling the B— Model" of Mariiio, Bouchard, Klemm, Pasquetti [8, 2] (at the moment for 
the refined topological vertex and the refined strip geometry). This seems to imply that mirror 
symmetry extends to the refined version of Gromov-Witten theory. 
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8. Examples 



8.1 The hexagon 

Consider the case where the 3 partitions bounding our plane partition are ah empty partitions 



A,/i,zv 



i.l) 



For the moment we consider that Smin = —a < 0, Smax = b > and A^ are finite. The profile of 
A = is: 

rx{s) = \s\/2, (8.2) 

and thus the domain for our non-intersecting paths is the hexagon bounded by the 6 lines: 



hexagon 



s = —a 
s = b 

h = -s/2 
h = s/2 

h = N + a + s/2 
^h = N + b-s/2 




(8.3) 



The corresponding matrix model of (4.15) reduces to 



Z 



dR (J'^Ri'^b-M-a) ^alntTrM-a ^{b-l)lnqTrMt ^^* ( *^ 'Voo ^^* ( ^ '^)oo 



(iHN) 



det (-t'^+i e-R; t)^ det {-q^ e^; q)^ 



(8.4) 

where M_a = diag(a/2, a/2 + 1, . . . ,a/2 + N - 1) and Mb = diag(V2, b/2 + 1, . . . ,b/2 + N - 1), 
i.e., 

b — a 



and 



Mb - M^a 



r^ ,. NiN - 1 Na 

Tr M_a = — ^ + 

2 2 



Id 



r^ ,. NiN - 1) iV6 

Tr Mb = — + 

2 2 



(8.5) 
(8.6) 



Our matrix model therefore reduces to a 1-matrix integral 

TrR det(-te-^;t)^ dei{-e^-q)^ 



Z (X 



h — a 

dR e^ 



{iHn) 



det (-t"+i e--^; t)^ det {-q^ e^; ( 



(8.7) 



The recipe of section 6.2 for finding the spectral curve, gives: 
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r t--(^) = XM = n{z) ^S^ 



-a< s<0 



-..(.) = (x,(z))i//^ = n(z)V/3 ,-/V/^y(.) , < . < 6 



1+Y{z) 



where ^(z) and y(z) must satisfy : 



u 
u 
u 
u 
u 
u 
u 



1 
1 

r 
r 



oo 
oo 



a-N 
b~N 



wheny = oo (8.8) 

whenl" = 

, whenl" = 

, wheny = oo 
wheny = — 1 
, wheny = —t"- 
, when y = —q 

The minimal assumption is to assume that the variable z lives in CP^, and u{z) and Y{z) have 
only the number of singularities required above: Y{z) has two poles and two zeroes, and u{z) has 
one pole and one zero of order 1, and one pole and one zero of order (3. Up to a reparametrization 
of z, we may assume that three of the points are 0, 1, oo, we make the Ansatz: 



Y{z) 



C 



{z-a){z--i) 
{z-5){z-e) 



u{z) = r 



1 



1 



w z 



p 



(8.9) 



The 7 unknown coefficients C,a,j,6,e,r,w are then determined by the conditions (8.8), which 

read in terms of this ansatz 

u{6) = 1 

u{j) = 1 

u{a) = t-^-" 

u{e) = t-^-^ (8.10) 

y(l) = -t" 

Y{l/w) = -q-^ 

y(o) = y(oo) = -1 

One can easily check that these 8 equations are not independent, only 7 of them are independent, 
and thus they allow us to determine the 7 unknowns. Plots of the corresponding arctic circle for 
various values of /3 generated by a Mathematica code are displayed in figure 13. 

8.1.1 Sending the size to oo 

The refined topological vertex is recovered in the infinite length limit of the hexagon. It turns 
out it is more convenient to start first to keep a and h finite and send iV — )■ oo. In this limit, 
Y{z) needs to have only one pole and one zero. Upon a reparametrization of z, we may choose 



The recipe for the spectral curve 



Y{z) = -z 

is rewritten as 

u{l) = 
it(0) = t-" 

u(po) = t"^ 

n(r) = oo 



^.11) 



^.12) 



u{q- 



oo 
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Figure 13: Numerical plots of the arctic circle for the hexagon with a ~ b ^ N, and [3 — 1,2, 3, 5, 10, 50. 
One side of the arctic circle is missing due to non real valuedness of the parameter in the parametric plot. 



that gives the spectral curve 




and we have 



XJz) 



Xsiz) 



f-b 



+-b-S- l-fz 



{l-zY+P 
(f^-z) {q->>-z)l^ 

l-q'z 
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s< 
s > 



^.13) 



^.14) 



See examples of plots of the limit shape in figure 14 for various values of a and b keeping (3 
fixed in the limit N ^ oo. Now it is straight forward to send a — t- oo and 6—7-00 and we end up 
with the limit of the spectral curve 



and we have 




{1-2)1+" 



Xsiz) 
Xs(z) 



t~2 



1-t^z 



t-l i^ 



zy 



s<0 
s> 



^.15) 



^.16) 
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Figure 14: with /3 = 2, we plotted the hexagons f" ^ t^ ^ 1/9, t" ^ £' ^ 1/100, t" ^ t^ ^ 1/900, 
f" =t^ = 1/1600 



This is the spectral curve of the refined topological vertex: 



spectral curve of the vertex 



y{z) = lny(2;) = Inz 

x{z) = h\u{z) = In - — — 



(8.17) 
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8.1.2 Other approach to the vertex's matrix model 

Now, let us choose a = b, and then send b — )■ oo, and write the matrix model of (8.4) 

2:0c / dR det{-te-^;t)^det{-e^;q)^ (8.18) 

which is the refined version of the matrix model obtained by Ooguri-Sulkowski-Yamazaki [31], 
also found by Sulkowski in [42]. 

Let us try to compute directly the spectral curve of the 1-matrix model (8.18). This should 
be easier, because it is a 1-matrix model, much better studied than the multi matrix model. How- 
ever, let us mention that we shall not get the same spectral curve, because our previous spectral 
curve was defined in terms of resolvents of matrices Mg, whereas here, we shall compute the 
resolvent of R. However, both spectral curves should be related by a symplectic transformation, 
and thus should have the same Fg^s, and thus compute the same partition function. 

The spectral curve of the general 1-matrix model j dRe~"^'^^^^\ is the large N limit of the 
resolvent 

/ 1 \^°^ 

M/(a;) =lnt(Tr -) (8.19) 

\ X — R I 

It is well known [32] that it lives on a double sheeted Riemann surface, i.e. 2 copies of the 
Riemann sphere glued along cuts, and along each cut it must satisfy: 

VF(x + iO) + l^(x-iO) = lnty'(a;), (8.20) 

Moreover, W{x) must be analytical in the physical sheet, and must behave like 

— I W{x)dx = N (8.21) 

•^''^ J support 

at large x in the physical sheet. 
Here we have: 

^(-+^o) + ^(--^o) = i-*Erv9^-i-^Erfl^ (8.22) 

k>\ k>0 ^ 

That gives in the small Int limit: 

W{x + iO) + W{x - iO) = In (1 + e^)'' (1 + e~^) + 0(ln t) (8.23) 

A solution is clearly 

^2W ^ {1 + Y)^+'' ^g_24) 

which is the same curve we have found in (8.17). Notice that in this case W{x) has only 
logarithmic cuts (it has no algebraic cuts as is customary when potentials are polynomial). 
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8.2 The cardioid 

Let us now consider A to be a rectangle, i.e., with 2 inner corners at vi and V2, and one outer 
corner m . We chose the origin of the s axis such that ui = ui + U2 , and we chose the origin of 
the h axis so that the upper corner of the bottom rectangle is at /i = (the bottom of the room 
is at a height h = 0). The domain is the following 





and thus the limit shape should have 3 tangents of slope +1/2, 3 tangents of slope —1/2 and 4 
vertical tangents. All the equations of section 6.2 can be easily written, and the spectral curve 
can be found at least numerically. In the unrefined case /3 = 1, the arctic circle is a degree 3 
algebraic curve, and when vi = —V2 = N/2, this is the famous regular cardioid often shown as 
an illustration of Okounkov's works [40] . 

However, let us further simplify by also chosing N = oo. Then, only the tangency conditions 
at the bottom of the domain remain, 




i.e., we need only one pole and one zero for Yi. We choose 



Yi{z) = -q^'-^' z , Y2{z) = -f^-^' 



z. 



(8.25) 



Then we have 



x.(z)=t"^x.,(.) ^-g,c::- 

Xs{z) = t^^ Xy2{z) ^Zt^2-^1^Z 

l-t''2-"i z 



Matching the 2 expressions at s = ui gives 



vi < s < Ul 

Ul < S < V2 
V2< S <U2 



(8.26) 



X,,{z) 



\-z 



\-q^ 



V'2 — fl 

t-^ X„,(z) 



\-z 



1 - t^a-"! Z 



^.27) 
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Then we want 



We thus find 




Xuo {z) 


= t-^ 


X,,{z) 


= t-^ 


Xui {z) 


= r^ 


Xv2 (z) 


1,2 


Xu^iz) 


U2 



X^,{q''^-'"^y/(^+'^^ =0 
X^,(t''i-«Og"i-^i) =00 
X,,(l)i//3 = oo 

Xy^il) = 00 



(l_g^i-«i^)/3 (l-t^2-«i2)/ 



(l-r2-«i^)/3 



(l-g^i-«iz)i+^ 



(1 - t'^o-viqvi-ui^-^ (1 _ _2)/3 (1 _ qU2~v2tv2-ui^y 



[l-q^i-^iz) (l-r2-«i2)/3 



(1 - t«0-t'lgt'l-«l^) (1 _ qU2-V2t'"2-Ui ^Y 

(l-t^2-«i^)l+/3 



(l-g^l^"lz) 



(1 - q«2-^2t^2-«l^)/3 (1 _ ^) (1 _ t«0-«l^fl-«l^) 



(1-z) (1 -t«0-flg^l-«l^) 

See a numerical plot of this spectral curve for /3 = 2 in figure 15. 

If then we send uq — )■ —00 and n2 — ?• 00, we find the spectral curve: 

Y {l + Yf 
This spectral curve is the should~be mirror of the geometry 




.28) 



(8.29) 



(8.30) 



8.3 Spectral curve for the general strip 

The above computation for the cardioid with one Ui and two fj's, can easily be extended to an 
arbitrary number of Uj's and Wj's. 

Again we work at A^ — )■ 00, so that Y needs to have only 1 pole and 1 zero, and we choose: 



Y,{z) 



^.31) 
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Figure 15: Plot of the limit shape obtained from the spectral curve (8.29) and with ji ~ 2 (and up to 
small artifacts from the numerics). 



and we recall that 
which we write 



y,+i(z) = r'+i-">g"»-^>y,(z), 



^.32) 



Yi{z) = -a z with Ci+i = T'+i-"' g"'-''' q. (8.33) 

similarly to the cardioid case, the minimal degree solution of all constraints of section 6.2 is: 



i.e., the mirror of the refined strip is: 



J J Cj z) 



iM) 



x = r2 



n,<. (1 + 1""^-'-"^ c, Y) n,>. (1 + 9"^-^^ c, yy 



with Q+1 = t'^'+l-"' g"«-^' Ci. 



^.35) 



9. Outlook &; Conclusions 



We construct the matrix model of the refined topological vertex and the refined strip geometry 
that is the building block of the superconformal SU{N) theory with Nf = 2N matter multiplets. 
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The corresponding spectral curves are found, thus the mirror of the corresponding geometries. 
Our matrix model differs from the previously proposed ones in the powers of the Vandermonde 
determinant, known as the /3— ensemble [43]. Unlike the previously constructed matrix models 
ours is an Hermitian one. The price to pay seems to be that we have a chain matrix model, and 
/? lies in the potentials instead of the Vandermonde. 

Having a matrix model, we first find its spectral curve, and we propose that this spectral 
curve should be the "mirror" of the corresponding refined geometry. Then, once the spectral 
curve is known, all string amplitudes are obtained by the topological recursion, in other words, 
the refined topological string also satisfies the "remodeling the B-model" construction. 

We find that the "mirror curve", i.e., the spectral curve of the matrix model, is a certain 
Legendre transform of the limit shape of the crystal model, in accordance with the unrefined 
case. For /3 G Q rational, the mirror curve is an algebraic curve. 

The strip geometry we consider in section 5. As mentioned two strips can be glued to give 
rise to the geometry that engineer superconformal gauge theories in Ad. In [12], using this fact 
and Gaiotto's construction of gauge theories combined with the AGTW relations the refined 
closed amplitude is connected to chiral A^v-i Toda three point function with one of the three 
primary fields of a special type [14]. Our approach, hence, give a matrix model interpretation of 
the three point function. 
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A. Appendix: Useful identities 

In this section we want to collect some of the identities and expressions we have made use of 
in the derivation of the crystal model for the refined strip geometry. When we have written 
the generating function for the plane partition we have encountered the following inner products 
which are written as shown in terms of Schur functions 



(7?|x^» nr_(l)t^»|a> = xHt^l'^ls./.Cl, t, . . . , e-'), (A.l) 

1=1 

a 

(M|x^«r+(1) n q'-'r+ilM = xl^^lg^l^U./^Cl, (?,..., g'^), (A.2) 

i=l 
a 

(Hx^«r_(i) nt'^''r_(i)|^) = xi-it^i^'is./^ii, t, . . .X), (A.3) 



j=i 



(7| n 9'^"r+(l)|i.) = q'^\^\s,/,{l, q,..., q^-'). (A.4) 



i=l 
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We have used the fohowing analytical continuation to express the Schur function with finite 
number of argument as Schur functions appearing in the topological vertex computations 

sa(1, t, . . . ; Qt~^, Qt~^ . . .) = sa(1, • • • , *"), with the choice Q = t"+^ (A.5) 

The Schur function on the left hand side of the above expression can be rewritten using 

s\{x, y) = ^ sx/n{x)s^{y), (A.6) 

where x and y stand for infinite many arguments. The following well-known identity has been 
used in the manipulations to show the equivalence 

where pi = —i + ^. 
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